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ABSTRACT 

Sco X-l has been the subject of many multi-wavelength studies in the past, being the 
brightest persistent extra-solar X-ray source ever observed. Here we revisit Sco X-l 
with simultaneous short cadence Kepler optical photometry and MAXI X-ray pho¬ 
tometry over a 78 day period, as well as optical spectroscopy obtained with HERMES. 
We find Sco X-l to be highly variable in all our datasets. The optical fluxes are clearly 
bimodal, implying the system can be found in two distinct optical states. These states 
are generally associated with the known flaring/normal branch X-ray states, although 
the flux distributions associated with these states overlap. Furthermore, we find that 
the optical power spectrum of Sco X-l differs substantially between optical luminosity 
states. Additionally we find rms-flux relations in both optical states, but only find 
a linear relation during periods of low optical luminosity. The full optical/X-ray dis¬ 
crete correlation function displays a broad ~ 12.5 hour optical lag. However during 
the normal branch phase the X-ray and optical fluxes are anti-correlated, whilst being 
correlated during the flaring branch. We also performed a Cepstrum analysis on the 
full Kepler light curve to determine the presence of any echoes within the optical light 
curve alone. We find significant echo signals, consistent with the optical lags found us¬ 
ing the discrete cross-correlation. We speculate that whilst some of the driving X-ray 
emission is reflected by the disk, some is absorbed and re-processed on the thermal 
timescale, giving rise to both the observed optical lags and optical echoes. 

Key words: accretion, accretion discs - X-rays: binaries - X-rays: individual: Sco 
X-l 


1 INTRODUCTION 

Scorpius X-l (hereafter Sco X-l) is a low-mass X-ray binary 
(LMXB), where a secondary main sequence star transfers 
material onto a magnetised neutron star (NS) via Roche- 
lobe overflow through an accretion disk. Sco X-l was the 
first, persistently bright, extrasolar X-ray source to be dis¬ 
covered (Giacconi et al. 1962). It lies at a distance of 2.8 


able nature across the electromagnetic spectrum (see e.g. 


kpc (Bradshaw et al. 1999), has an orbital period of 18.9 


hours (Gottlieb et al. 1975 Cowley & Crampton 1975 Hynes 


& Britt 2012), and contains a O.4M0 M-type companion 
(Steeghs & Casares 2002). Since its discovery, many studies 
have been pursued on Sco X-l, all revealing its highly vari- 
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Mook et ah1|1975| |Hertz et al.||1992| |Dieters &; van der Klis| 

2000| ). In particular, Sco X-l is a member of a group of 
objects called Z sources. The class name arises due to the 
fact that apparent Z-shaped tracks appear when soft X-ray 
colours are plotted against hard X-ray colours. Many mul¬ 
tiwavelength studies have been pursued on Z-sources in the 
past, revealing how the mass accretion rate, ultraviolet line 
and continuum fluxes and optical brightness are all corre¬ 


lated with where the sources lie on the Z-track ( Hertz et al. 
1992 Hasinger et al. 1990 Vrtilek et al. 1990[ 1991). The 


top-most portion of the Z-track is referred to as the horizon¬ 
tal branch (HB), and systems have been shown to be rarely 
found in this state. As the sources progress downward, they 
enter the normal branch (NB), whilst during episodes of 
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higher X-ray flux, they enter the bottom part of the Z-track 
called the flaring branch (FB). In the case of Sco X-l, the 
HB is far less pronounced than either the NB or FB. Ad¬ 
ditionally to X-ray colour changes the X-ray power spectral 
density (PSD) of Sco X-l displays clear quasi-periodic os¬ 
cillations (QPOs), which appear to increase in frequency as 


the source moves downwards along the NB ( |van der Klis 
et al.||1996| . 


One of the first models invoked to explain the Z-track 
phenomenology was that of Priedhorsky et al. (1986) and 


Psaltis et al. (1995). In that model most of the X-ray emis¬ 


sion originates from the NS magnetosphere, the hot central 
corona, and an extended corona through which material falls 
radially towards the NS. Variations within the Z-track then 
reflect variations in the mass accretion rate which varies be¬ 
tween 0.5—1.1 times the Eddington critical rate, Me- During 
low accretion rate episodes, Sco X-l is found in the NB and 
moves downwards as M increases. This causes increased ra¬ 
diation pressure allowing material to pile-up close to the NS. 
The optical depth then increases, causing more X-ray pho¬ 
tons to be absorbed and re-emitted as optical photons. This 
phase is associated with an X-ray to optical anti-correlation 
( McNamara et ah] 2003). As the source transitions between 
the NB and the FB, the mass accretion rate exceeds the crit¬ 
ical Eddington rate, causing the mass flow to become chaotic 
with material in the disk flowing both radially inwards and 
outwards. As material keeps building up in the outer corona 
the electron scattering opacity increases, redirecting some of 
the X-ray emitted photons onto the accretion disk. Contrary 
to the NB phase, the FB phase produces correlated X-ray 
and optical fluxes. 

Recently, a different model has been proposed based on 


the idea of an extended accretion disc corona (ADC, Church 


V Balucinska-Church 1995; Church & Baluciriska-Church 


2004 


Church et al.|2012 ). In this model M increases as the 


source moves upward along the NB track (the opposite of the 
Psaltis et al.| l 995| model). Although qualitatively different, 


this would still predict an X-ray/optical anti-correlation due 
to increased optical depths with higher accretion rates. Fur¬ 
thermore, movements along the FB away from the NB/FB 
vertex are also associated with increased M, but here the 
flaring behaviour is associated to unstable nuclear burning 
on the NS surface (Church et al. 2012). The M increase 
along the FB is, however, only a property of Sco X-l like 
sources: other systems like Cyg X-2 are instead thought to 
have constant M along the FB. Thus, given that most of 
the X-ray radiation is released from the NS rather than the 
disc, we might also expect that this scenario would produce 
correlated X-ray and optical fluxes as in the previous model. 
Furthermore, as movements along both FB and NB branches 
are associated with increased mass-transfer rates (and thus 
increased radiation pressure), the model predicts that both 
should display jet-launching. Indeed both branches have 
been shown to display radio emission, although the NB is 
more radio-loud than the FB ( Hjellming et al.|1990 ). 


Bradshaw et al. (2007) have studied the correlations 


between the X-ray spectral characteristics and QPOs in Sco 
X-l. They found a strong correlation between the kHz QPOs 
and the spectral power law index and interpret this as be¬ 
ing due to changes in the geometrical configuration of the 
corona. The same study shows how the equivalent width 
(EW) of the Iron K a emission line is anti-correlated with 


spectral hardness. The increase in EW implies a stronger 
accretion disk wind, assuming the Iron K a line is generated 
in the wind itself. 

With the aim of testing and refining the current mod¬ 
els for Sco X-l, we revisit this system with simultaneous 
data obtained from the K2 mission and the MAXI instru¬ 
ment, augmented by ground based optical spectroscopy ob¬ 
tained with HERMES. Our data reduction procedures are 
described in Section [2] whilst our results are presented in 
Section [3] Our discussion, together with interpretation of 
our results, is discussed in Section [4] whilst our conclusions 
are drawn in Section [5] 


2 DATA REDUCTION AND ANALYSIS 

In this section we present the three datasets used in this 
work (K2, MAXI , HERMES ), and describe any relevant 
data reduction procedures employed. 


2.1 K2 light curve 


The NASA Kepler mission was successfully launched on 
March 7 2009. It continuously pointed towards the same 116 
square degree field-of-view (FOV) during the entire mission, 
and obtained light curves with a 58.8 second cadence (short 
cadence, SC) or with a 30 minute cadence (long cadence, 


LC) during every quarter (Q) of observations (Boruck i et al.| 
2010). In May 2013, Kepler entered a prolonged safe mode 


due to a second reaction wheel failure. Since then the mis¬ 
sion has been re-purposed as the K2 mission, covering new 
areas of the sky with the same SC and LC modes and the 
same FOV size. K2 thus stares for 3 months at the same 
field before moving onto the next one. Although K2 is not 
able to obtain year-long light curves as the original Kepler 
mission, it does allow to study the timing properties of a 
wider range of astrophysical objects as the mission sweeps 
through the so-called Kepler ecliptic plane. 

Sco X-l was observed by the K2 mission during Cam¬ 
paign 2 between August 23 and November 13 2014 (78.8 
days), and is the first ever LMXB to be observed with Ke¬ 
pler. Sco X-l has a registered Kepler magnitude (Kp) of 
12.4 in the K2 Ecliptic Plane Input Catalog (EPIC). Here 
we analyse SC data (cadence of 58.8 seconds) obtained from 
the Mikulski Archive for Space Telescope (MAST) archiv^] 
The data is provided in raw format, consisting of target pixel 
data. For each 58.8 second exposure we thus have a 14 x 
12 pixel image centered on the target. Sco X-l was near the 
edge of module 15.3. Although this module is not known to 
be affected by Moire Pattern Drift (MPD) noise the target 
point spread function (PSF) is asymmetric as it lies at the 
edge of one of the outermost modules. As no other sources 
were present within the 14 x 12 pixel images, we created the 
light curve by manually defining a large target mask as well 
as a background mask. A large target mask is required due 
to occasional small scale jittering of the spacecraft, result¬ 
ing in the target moving slightly from its nominal position. 
Fig. [l] shows the average target pixel image obtained from 


1 http://http://archive.stsci.edu/k2/ 
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Figure 1 . K2 average target pixel image for Sco X-l. The target 
mask pixels are marked in red, whilst background mask pixels are 
marked in black. 
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Figure 4. Sco X-l colour-colour X-ray diagram obtained from 
MAXI data. The soft colour is defined as the ratio of 4-10 keV 
/ 2-4 keV, whilst the hard colour is defined as the ratio of 10-20 
keV / 4-10 keV. Typical error bar shown on the top-left. 
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2.2 MAXI light curve 


Since the MAXI (Monitor of All-sky X-ray Image; Matsuoka 
et al.||2009| experiment on board the International Space 
Station (ISS) started in 2009 August, the GSC (Gas Slit 
Camera; Mih ara et al.|20lT Sugizaki et al.|20TT ), one of the 
two MAXI detectors, has been scanning almost the whole 
sky every 92-minute orbital cycle in the 2-30 keV band. Here 
we use the MAXI light curve of Sco X-l during the period 
when the system was also observed by the K2 mission. We 
obtained the light curve through the MAXI websit^J which 
provides reduced light curves in the three energy ranges 2- 
4 keV, 4-10 keV and 10-20 keV. The combined light curve 
is shown in the bottom panel of Fig. [2] The correspond¬ 
ing X-ray colour-colour diagram is shown in Fig. [4] clearly 
displaying the NB/FB dichotomy. 


2.3 Mercator/HERMES spectra 


Figure 3. Fourier Transform of the full 78.8 days SC K2 data 
for Sco X-l, showing a clear peak at the orbital period of 18.9 
hours. 


115,679 individual target images. We removed 1,919 obser¬ 
vations because of bad quality due to occasional spacecraft 
rolls or due to cosmic rays. Fig. [I] also shows in red the 
target mask and in black the background mask. We pro¬ 
duce the lightcurve by summing together all target pixels 
for each exposure, and subtract the average background ob¬ 
tained from the background pixel mask. The obtained light 
curve is shown in the top panel of Fig. [2] We believe this 
is the best ever optical light curve obtained for an LMXB 
in terms of timespan, cadence and photometric quality. A 
Fourier transform of the whole light curve reveals the or¬ 
bital period at 18.9 hours (see Fig. |3|. 


We observed Sco X-l with the HERMES spectrograph 
(Raskin et al.||2011 ), mounted on the 1.2 m Mercator Tele¬ 
scope at La Palma, Canary Islands, Spain. This highly ef¬ 
ficient echelle spectrograph has a resolving power of R = 
86000 over the range 3800A to 9000A. Each spectrum was 
obtained from a 45 minute exposure. The raw spectra were 
reduced with the instrument-specific pipeline, but is not 
flux-calibrated. In total we obtained 6 spectra at the times 
marked with the vertical dashed lines in Fig. [2] Here we 
specifically examine the evolution of the 6 emission lines Ha 
6563 A, H 0 4861 A, Hel 5876 A, Hel 6678 A, Hel 7065 A 


and Hell 4685 A. To do this we normalise each line with an 
estimate of the neighbouring continuum. For each line, we 
iteratively fit a first degree polynomial around 40A of the 
desired line centre using sigma-clipping. We found this to 
converge after about 5 iterations. In cases where no evident 


2 http://maxi.riken.jp 
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Figure 2. Top panel: K2 Sco X-l light curve. The system was observed for over 78 days at 58.8 second cadence. The units on the y-axis 
are electrons/second, and can be converted to Kepler magnitudes Kp using the conversion found in the Kepler Instrument Handbook. 
Bottom panel: 2-20 keV MAXI light curve during the same period. Data points are colour coded using the NB/FB decomposition shown 
in Fig. [T] Typical MAXI error bar is shown on the top-right. The vertical dashed lines mark the times of our HERMES observations. 


emission line was present, the iterative sigma-clipping pro¬ 
cedure excluded ~ 3% of data, whilst in cases where clear 
emission lines were present it excluded ~ 15%. The resulting 
normalised line profiles are shown in Fig. [5] together with 
their measured equivalent widths. 


3 RESULTS 

In this Section we investigate the properties of Sco X-l in 
the context of its NB and FB track dichotomy. We will de¬ 
termine what method best separates the 2 tracks using the 
K2 and MAXI datasets and present average PSDs for the 
low and high optical luminosity states. This section addi¬ 
tionally presents the results of cross-correlating the MAXI 
and K2 light curves, the Cepstrum analysis to detect pos¬ 
sible echoes in the K2 light curve, as well as our optical 
spectroscopy campaign results. 

3.1 Normal and Flaring branch properties 

Fig. [6] shows the optical flux distribution obtained from the 
K2 light curve (left panel, gray line) and the X-ray flux dis¬ 
tribution obtained from the MAXI light curve (right panel, 
gray line). The K2 flux distribution is clearly bimodal, and 
possibly the MAXI fluxes as well. From the flux distribu¬ 
tions alone it is not clear whether the bi-modality is caused 
by the NB/FB dichotomy. To determine whether this is the 
case we performed principal component analysis (PCA) de¬ 
composition on the MAXI X-ray colours shown in Fig. [4] 


in order to determine which observations belong to which 
Z-track branch. Fig. [7] shows the obtained PCA projection. 
By using principal component 2 (PC 2) we can clearly sepa¬ 
rate the NB from the FB. We employ conservative ranges to 
exclude transition data points and select FB MAXI obser¬ 
vations as having PC 2 >0.2 and NB MAXI observations as 
having PC 2 < —0.2. The resulting X-ray flux distributions 
for both the NB and FB obtained from our selected points 
are shown in Fig. [6] Although NB X-ray fluxes seem to only 
populate the low flux distribution in Fig. [6] the FB fluxes 
seem to populate both distributions. Thus X-ray fluxes alone 
are not a good discriminator between NB/FB tracks. We 
then located all K2 observations that lie within ±45 min¬ 
utes from the MAXI observations (MAXI data points have 
a cadence of one ISS orbit, or 1.5 hours). The resulting FB 
and NB flux decomposition from our PCA projection is also 
shown in Fig. [6] It is clear that in the absence of X-ray colour 
information, optical fluxes are better than X-ray fluxes in 
distinguishing between NB and FB. However we point out 
that optical fluxes alone are also not good indicators of the 
NB/FB dichotomy, given the long tails displayed in Fig. [6] 

Given the optical flux bimodality, we decided to pro¬ 
duce broad-band PSDs from the K2 light curve for differ¬ 
ent flux levels. We selected conservative ranges to avoid the 
transition region, and chose to call low optical luminosity 
observations all data lying below 4.6 x 10 4 e _ s -1 and high 
optical luminosity those lying above 5.6 x 10 4 e _ s -1 . These 
limits are marked with vertical lines in Fig. [8] We then com¬ 
puted 1-day, non-overlapping PSDs (78 of them), and pro¬ 
duced a low-luminosity PSD by averaging 30 of them whose 
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Figure 5. HERMES spectra normalised to continuum. Each row is one observation taken at times marked in Fig[2]with vertical dashed 
red lines. Columns are different zooms around different emission lines. Each panel also shows the measured equivalent width of the lines. 
Note the increased line equivalent width when Sco X-l was in its Normal branch state (Spectra 3 & 4). 




Figure 6. Left panel: K2 flux distribution decomposed using PCA as applied to X-ray colours. Note that due to the conservative range 
employed the Flaring and Normal branch, data points are not all recovered. Right panel: MAXI flux distribution decomposed using PCA 
as applied to the corresponding X-ray colours. 
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Figure 7. Projected Hard and Soft X-ray colours using principal 
component analysis. We define conservative regions for the Flar¬ 
ing and Normal branch by selecting MAXI observations that lie 
0.2 above/below principal component 2. 



Figure 8. Sco X-l optical flux distribution obtained with K2. 
The vertical lines correspond to conservative boundaries for low 
and high optical luminosity states determined by eye. 


mean fluxes lay within the optical flux range shown in Fig. 
[8] The high optical luminosity PSD was similarly obtained 
by averaging 28 PSDs. The logarithmically binned result of 
this procedure is shown in Fig. [9] At first look it is already 
clear that the PSD shapes differ substantially between the 
low and high optical luminosity states. We include in Fig. [9] 
a dashed gray line showing a power-law with a —2 gradient 
for comparison. During both high and low optical luminos¬ 
ity states the PSDs appear to follow a general —2 power-law 
shape. At low optical luminosities the PSD possibly displays 
an additional noise component around ~ 2 x 10 -4 Hz. On 
the other hand during high optical luminosities a clear broad 
noise component is observed in addition to the general ~ — 2 
power-law slope, peaking around ~ 10 _3 Hz. The flattening 
observed at the highest frequencies is probably due to in¬ 
strumental noise. It is clear that the PSD phenomenology 
is quite complex and differs substantially between states. A 
detailed modelling in terms of power laws and Lorentzian 
components is beyond the scope of this work and will be 
addressed in future. 

Using the optical flux limits defined in Fig. [8] we are 
also able to produce root-mean-square vs. flux (rms-flux) re¬ 
lations for both the low and high optical luminosity states. 
This was done by computing mean fluxes and rms values 
for all non-overlapping 5-minute data segments, excluding 
any data gaps caused by the flux limit selections. Fig. [To| 
shows our binned results. In the low-luminosity state, an 
rms-flux relation consistent with linearity is obtained. In 
the high-luminosity state, it appears that linearity is bro¬ 
ken at the highest fluxes. This could be suggestive that an 
additional flux distribution exists at the highest fluxes, thus 
making the light curve non-stationary and breaking the rms- 
flux linearity. We note that this could be connected to the 
trimodal flux distribution previously reported in Sco X-l 
( Mook et al.|1975 McNamara et al.|2003| and/or the extra 
variability component observed during the high optical lu¬ 
minosity state as seen in Fig. [9] A similar property might 



Figure 9. PSDs obtained from the low and high optical lumi¬ 
nosity intervals defined from the flux distribution shown in Fig[8] 
The dashed gray line shows a —2 power-law for comparison. 


also be observed in the NS XRB PSR J1023+0038 (Bog¬ 


danov et al. 2014), although we caution the reader that 


this particular system might have a higher accretion rate 
that Sco X-l, and thus any phenomenological similarity in 
the flux distributions might be coincidental. Nevertheless, it 
would be interesting to establish whether the rms-flux rela¬ 
tion holds its linearity for this system or whether it breaks 
down at high fluxes as seen in Sco X-l. We also note that at 
the lowest fluxes the high optical luminosity rms-flux rela¬ 
tion gradient is higher compared to the low luminosity one. 

Fig. [5] shows the 6 spectra obtained with the HERMES 
spectrograph. Each column shows 6 different continuum nor¬ 
malised emission lines. The low ionisation emission lines of 
Hcq H f3 and Hel qualitatively appear to have larger equiva- 
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Figure 10. rms-flux relations obtained for Sco X-l during low and high optical luminosity states selected using the regions defined in 

Fig m 


lent widths during periods of low optical luminosity (Spec¬ 
tra 3 & 4) compared to the high luminosity spectra. On the 
other, hand the equivalent width of the higher ionisation 
Hell line does not appear to change by much. Furthermore, 
Hell appears to be double-peaked only during the high op¬ 
tical luminosity states phase. Possible explanations for this 
qualitative behaviour will be discussed in Section^] 


3.2 X-ray lags &; optical echoes 

We computed the discrete cross-correlation function (DCF, 
Edelson fc Krolik||1988 ) between the K2 and MAXI 2-20 
keV lightcurves using a time step of 45 minutes (half the 
MAXI cadence). The results are shown in Fig. 


11 


where the 

x-axis is optical lag in units of hours. We further computed 
confidence levels by simulating 10 5 lightcurves on the same 
K2 sampling pattern of the original data. Each light curve 


was simulated using the Timmer X Koenig (1995) algorithm. 


As we do not have an analytical model for the intrinsic PSD 
of Sco X-l within the frequency range of interest, we pro¬ 
duced the input PSD by first computing the full Fourier 
transform of the original K2 data and smoothing it with a 
50-point moving average. The dashed line in Fig. EH shows 
the obtained 3cr (99.7%) contour level from our Monte-Carlo 
simulation. 

A clear significant zero-lag peak is evident. Given our 
timing resolution, this peak is consistent with the find¬ 
ing that the optical light curve lags the X-ray light curve 
by tens of seconds due to reprocessing from the sec¬ 


ondary star and/or accretion disk (Munoz-Darias et al. 
2007 Britt| 2013). Additionally the DCF is very asymmet¬ 
ric, with a second broad significant bump around 12.5 hours. 
This timescale might be consistent with the disk thermal 
timescale, and will be discussed further in Section [4] 

Motivated by the different properties between the NB 
and the FB, as well as the idea that X-ray/optical light 
curves should be anti-correlated in the NB, we further com¬ 


puted the DCF for the two separate populations. We used 
the full K2 lightcurve for both, and selected MAXI data 
segments using the PCA decomposition shown in Fig. [7] 
The results of this are shown in Fig. |12| where the contour 
levels have been computed in the same way as for the whole 
dataset of Fig. EH Although slightly more scattered, the 
DCFs for the individual NB and FB show different quali¬ 
tative behaviour. Similarly to the full dataset DCF, the FB 
DCF is asymmetric in the same direction. It displays a pos¬ 
itive peak at zero time-lag, and although marginally signifi¬ 
cant, a potential secondary optical lag peak at ~ 4.5 hours. 
During the NB the two light curves are anti-correlated and 
appear more symmetric than for the FB. This effect can 
also be noticed in the light curves themselves by close visual 
inspection (see e.g day « 75 in Fig.[2|. The cor relation/anti¬ 
correlation properties between X-ray and optical fluxes have 
been already reported in the past (McNam ara et al.| |2003) 
using X-ray flux vs. optical flux plots, which we here con¬ 
firm and characterise in more detail using the DCF. Fig. |13| 
shows our version of X-ray vs. optical flux diagram, where 
we have additionally marked with circles the positions of 
when our optical spectra was taken. Simply from visual in¬ 
spection it can be noticed that Spectrum 3 & 4 (which have 
the strongest emission lines) have been taken at times where 
the optical flux was lowest, but where the X-ray flux is not 
so different as for the other spectra. This suggests that low 
ionisation line strength is more strongly correlated to optical 
flux and not X-ray fluxes or colours. 

Finally we were also interested in determining whether 
the K2 optical lightcurve alone displays any echoes. Any de¬ 
tection of such features would be interesting for a number of 
reasons. For example, if the optical lag (response) detected 
in Fig. is associated with the disk thermal timescale, we 
might expect an echo within the optical light curve alone 
around the same time. To determine whether any echoes 
are present in the K2 light curve we computed the so-called 
Cepstrum. The Cepstrum is computed as the inverse Fourier 
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Normal branch Flaring branch 




Optical lag (hours) Optical lag (hours) 


Figure 12. Discrete cross-correlation between MAXI and K2 data selected according to the PC A decomposition for normal and flaring 
branches (solid line, left and right respectively). The dashed lines mark the 3cr detection levels. 
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Figure 13. K2 vs. MAXI fluxes colour coded using the NB/FB PCA decomposition of Fig. [t] Circles indicate the time when the 6 
HERMES spectra were obtained, numbered from 1 to 6 as in Fig. [5] 
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Figure 11. Discrete cross-correlation between MAXI and K2 
data (solid line) and 3a detection level (dashed line). 


transform of the logarithm of the estimated spectrum of the 
light curve (see e.g. Bogert et al. 1963| Afraimovich 1981 


Maccarone et al.|2015 ), and is used to find repeating signals 


(echoes) within a single time-series. Confidence limits can 
also be computed with a similar method as that employed 
for the DCF confidence limits. We simulate 10 5 light curves 
using the same input PSD as previously adopted, and de¬ 
termine for each timescale the 3cr level. The results of this 
is shown in Fig.[l4] Three significant signals are detected at 

1 hour, 4.4 hours and 12.6 hours. We will discuss possible 

mechanisms to generate these echoes in Section [4] but al¬ 
ready note here that the signal at 12.6 hours is consistent 
with the observed optical lag shown in Fig. eld and that the 
4.4 hour signal is consistent with the FB DCF marginally 

detected secondary peak in Fig. [12] We also point out that 
the Cepstrum is solely based on the K2 light curve, thus 
any matching signals between the Cepstrum and the DCFs 
are unlikely to be coincidental. 


4 DISCUSSION 


In this Section we discuss our findings in the context of the 


Church et al. (2012) model for Z-sources in general and Sco 


X-l in particular. We also propose interpretations for several 
newly discovered phenomena. 


4.1 X-ray and optical state properties 


Phenomenologically, Sco X-l transitions between NB and 
FB on timescales of « 1 day. According to |Church et al.| 


(2012) the proposed mechanisms for generating these tran¬ 


sitions is related to the mass-transfer rate varying, such that 
M increases along both the FB and NB tracks away from the 
vertex. The X-ray colour changes reflect the relative contri¬ 
bution to the X-ray flux of the black-body disc component 
(dominating in the soft X-ray band) and comptonised corona 
(dominating in the hard X-ray band). As the optical emis¬ 
sion is thought to originate from a colder outer disk, the op¬ 
tical flux distribution better tracks M changes, being unaf¬ 
fected by X-ray absorption and/or the X-ray contribution by 
different X-ray emitting components. Thus, although both 


the X-ray and optical flux distributions appear bimodal, the 
bimodality found in X-rays may not correspond simply to 
states with different M. 


One important consequence of the 

Church et al. (2012) 

model (but also the older Psaltis et al. 

1995 model) is that 


as M increases, during the NB phase, the optical depth close 
to the NS also rises, due to higher radiation pressure. This 
might cause X-ray photons to be absorbed and re-emitted 
as optical photons, resulting in an anti-correlation between 
X-ray and optical fluxes. This is observed in Fig. m On the 
other hand, during the FB phase, the Church et al. (2012) 


model predicts an additional X-ray component caused by 
unstable nuclear burning on the NS surface. This additional 
X-ray component would shine on the optical emitting outer 
accretion disk, so we should observe a correlation between 
optical vs. X-ray fluxes. This is also observed as shown in 

Fig. [12] 

We have also found that the broad-band PSD shape 
of Sco X-l differs substantially between the low and high 
optical luminosity states, with the latter displaying an addi¬ 
tional broad noise component breaking at : 


10 3 Hz super¬ 


imposed ~ —2 power-law. Bildsten (1993) has previously 


reported that a possible time-dependent accretion regime 
exists in accreting NS, where small patches of the NS burn 
intermittently (unstable and persistent, as opposed to burn¬ 
ing steadily or rapidly on the entire NS surface). The fire 
will then propagate around the NS on the timescale it takes 
to accrete enough fuel for the next instability, which lies in 
the range between 10 * 1 * 3 - 10 4 seconds. Thus, low-level lumi¬ 
nosity variations should be observed on a similar frequency 


range. Indeed Hasinger & van der Klis (1989) have reported 


“very low frequency” noise observed on similar timescales. 
Given that the high optical luminosity state is somewhat 
associated to the FB for Sco X-l, and given that accord¬ 


ing to the Church et al. (2012) model this is where nuclear 


burning occurs on the NS surface, it is possible that the addi¬ 
tional noise component in the high optical luminosity state 


observed here is in fact associated to the Bildsten (1993) 
mechanism. 

It is also interesting to note how the rms-hux rela¬ 
tions differ between the two low and high optical luminos¬ 
ity states. Linear rms-hux relations have been detected in 
ah accreting compact objects (from Active Galactic Nu¬ 
clei to accreting white dwarfs, see e. g |Uttle y et al.| 2005| 
Heil et aT7||2012| |Scaringi et al.||2012[ |Van de Sande et al.| 
2015| ). The general model invoked to explain the linear rms- 
hux relations is that of the huctuating accretion disk model 


( Lyubarskii|1997[ Kotov et al.|2001 Arevalo &; Uttley|2006 ) 
In this model, the observed variability is associated with 
modulations in the effective viscosity of the accretion how 
at different radii. More specifically, the model assumes that, 
at each radius, the viscosity - and hence accretion rate - 
huctuates on the local viscous timescale around the mean 
accretion rate, whose value is set by what is passed inwards 
(again on the viscous timescale) from larger radii. The over¬ 
all variability of the accretion rate observed through light 
curve variations is therefore effectively the product of ah 
the huctuations produced at larger radii. It is important to 
realise, however, that the detection of the rms-hux relation 
only suggests that the observed variability is the result of 
multiplicative processes (rather than additive) and does not 
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Figure 14. Cepstrum obtained from the K2 light curve (solid line) and 3a contour level (dashed line). The x-axis of the Cepstrum is 
referred to as the “quefrency”, and has unitis of time. Three significant peaks are detected at 1 hour, 4.4 hours and 12.6 hours. 


necessarily rule out more complex (or less complex) accre¬ 
tion disk models to generate the observed variability. 

In Sco X-l, the optical K2 light curve does always 
exhibit a positive correlation between mean flux and rms. 
However, this correlation is only linear in the low optical 
luminosity state. This suggests that during low optical lu¬ 
minosities the light curve is indeed the result of multiplica¬ 
tive processes, which are possibly the result of a fluctuating 
accretion disk. During times of high optical luminosity the 
rms-flux relation deviates from linearity (with an apparent 
change of slope to lower values) suggesting that more com¬ 
plex phenomena are at play during this phase. It is possible 
that the high optical luminosity state is composed of more 
than one accretion regime (for example the unstable nuclear 
burning mechanism of Bildsten 1993), causing the higher 


fluxes of this phase to be non-stationary with respect to the 
lower fluxes. Whatever the reason, it is interesting to note 
that multiplicative processes alone cannot reproduce non¬ 
linear rms-flux relations. 


Fen der et al.| 2009). In this scenario the optical continuum 
would drop much faster than the ionizing continuum pro¬ 
duced in the inner disc regions. The EW of the optical lines 
produced in the outer disc should thus increase. The same 
would also hold during the opposite transition from low to 
high states, as the continuum would rise much faster than 
the emission lines, thus swamping them out and decreasing 
their measured EWs. Another possibility is that during the 
high optical luminosities most of the accretion disk becomes 
too hot for the production of low ionisation emission lines. 
This would not necessarily cause higher ionisation lines such 
as Hell to decrease in EW. In this case optical fluxes alone 
(rather than X-ray colours) would provide a better tracer for 
M changes in the optical emitting region. Indeed 2 of our 
6 spectra with the lowest optical fluxes display the highest 
EWs. Furthermore, we do not measure significant changes 
in the Hell EW, possibly making this scenario a more viable 
explanation. 


4.2 Optical spectroscopy 

The HERMES spectra reveal that the low ionisation lines 
of Ho, H/3 and Hel appear stronger when Sco X-l is in a 
low optical luminosity state (but not necessarily at low X- 
ray luminosities). |Fen d er~3t al.| (2009) have found a simi¬ 
lar trend in a sample of XRBs, where an anti-correlation 
is found between the X-ray luminosity and the Ho equiva¬ 


lent width. Fender et al. (2009) suggest that the observed 
anti-correlation is the result of optical-depth changes. 

In Sco X-l, the anti-correlation could be due to the 
outer disk becoming cold and optically thin in the con¬ 
tinuum, but not the emission lines, during the transition 


Although the Hell EW does not appear to significantly 
change between the low and high optical luminosity states, 
it does display a double-peaked profile only during low lu¬ 
minosity ones. We note that the Hell line profile changes do 
not appear to depend on the NB/FB dichotomy, but simply 
on optical luminosity (see Fig. [Tdj ). The change from dou¬ 
ble to single peaked profiles has been previously associated 
with the presence of outflows (see e.g. Murray & Chiang 
1996). This deduction of a wind could potentially be con¬ 


between the two optical luminosity states (Williams 1980 


sistent with the X-ray observations of the Iron K a emission 
line modelling in Sco X-l (Bradshaw et al.| 2007), which as¬ 
sociates a strong wind to an X-ray softening. 
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4.3 X-ray time-lags and optical echoes 


One other interesting feature is the broad ~ 12.5 hour opti¬ 
cal lag shown in Fig. pT] Naively, from the fluctuating accre¬ 
tion disk model (Lyubarskii 1997 Kotov et al.|2001 Arevalo 


fe Uttley| 2006 ), we might have expected to observe an X- 


ray lag due to material propagating from an optical emitting 
region in the outer disk to an X-ray emitting region in the 
inner disk. Instead we observe the exact opposite for Sco 
X-l: first we observe X-rays and hours later the signal ap¬ 
pears in the optical band. Given the hour-long measured 
timescales in the cross-correlation function of Fig. mi our 
result rules out light-travel time reverberation from the disk 
as a possible mechanism. We do not know at the time of 
writing what is the mechanism behind the observed optical 
lags but speculate that it could be the result of reprocessed 
X-ray photons to optical wavelengths on the disk thermal 
timescale. This idea is not new, and has been invoked to 
explain Fourier-dependent time-lags observed in accreting 
white dwarfs (Scaringi et al. 2013). 

Although the optical emission in the broad Kepler pass- 
band will originate from a large radial disk component, we 
can attempt to place constraints on a since the thermal 
timescale tth = ^ (where is the dynamical frequency at 
some given radius in the disk and a is the viscosity parame¬ 
ter of Sh akura k, Sunyaev|19 73). Fig. 15 shows the obtained 
a values as a function of disk radius assuming a NS mass of 
1.6M 0 ( |Munoz-Darias et al.| |2007) by setting t t h — 12.5 or 
Uh = 4.4 hours (two of the peaks detected in the Cepstrum 
analysis). The constraints obtained from Fig. 15 are physi¬ 
cally plausible, suggesting that the measured lag timescale 
could be associated to the disk thermal timescale. 


If the measured lags are really associated to the thermal 
timescale we would expect the X-ray and optical variability 
to have a very small, light-travel time, lag which we ob¬ 
serve in Fig. m as the strong zero-lag peak. As some X-ray 
radiation will be absorbed by the disk and re-emitted on 
the thermal timescale at optical wavelengths we might then 
also expect a secondary, broad, optical lag peak. In this case 
we might expect the optical lightcurve to also display an 
echo, since the same signal should appear twice within the 
lightcurve (once reflected and once reprocessed), with a de¬ 
lay close to tth • As one of the significant peaks found in the 
Cepstrum (Fig. p~4| is similar to the measured broad X-ray 
lag, we tentatively associate this with the echo mechanism 
proposed above. A similar process might explain the other 
2 detected echoes. In fact, the FB DCF in Fig. [T2] is quali¬ 
tatively similar to the full DCF of Fig. [Tl] Given our time 
reolution used to compute the DCFs, we would not be able 
to observe the 1 hour cepstral peak in the DCFs. We point 
out that our Cepstrum analysis is based on the full 78-day 
K2 light curve, and that we are not able at this time to 
localise the echoes within the light curve. Thus multiple fea¬ 
tures in the Cepstrum do not necessarily need to occur at 
the same time within the 78-day period. 

We caution at this stage that our tentative association 
of the optical lags and Cepstrum signals to the disk thermal 
timescale are tentative, as they have never been observed 
on these timescales. However, given the long timescales in¬ 
volved, the non-detection of such lags in other similar sys¬ 
tems is not surprising, since one would require continuous 
and long monitoring of XRBs at both X-ray and optical 



Figure 15. a dependence with disk radius for t t h ~ 12.5 hours 
(solid line) and tth ~ 4.4 hours (dashed line). We also mark ap¬ 
proximate radii regions of different electromagnetic emission. 


wavelengths. It is interesting to note the possibility that op¬ 
tical lags might only be detected in the FB state (since the 
NB shows X-ray/optical anti-correlation), and that the disk 
heating mechanism could be driven by the unstable nuclear 
burning on the NS surface. A detailed study into whether 
this scenario is feasible for Sco X-l is however beyond the 
scope of this paper. 


5 CONCLUSION 

A prime motivation for this study has been to revisit the 
well-known LMXB Sco X-l with newly available data to 
test and refine the current physical model used to explain 
its phenomenological behaviour. Our dataset comprises a 
continuous 58.8 second cadence light curve over more than 
78 days obtained with the K2 mission, and represents the 
best ever optical light curve obtained for a LMXB. We ad¬ 
ditionally used X-ray data obtained with the MAXI instru¬ 
ment on-board the ISS which provided us with X-ray light 
curves in the three energy bands 2-4 keV, 4-10 keV and 10- 
20 keV over the same observing period of K2. Finally, we 
augmented our analysis with high resolution optical spec¬ 
troscopy taken during the same period with the HERMES 
instrument mounted on the 1.2 meter Mercator telescope. 

We can summarise the major findings of this paper as 
follows: 


(i) In the absence of X-ray colours, optical fluxes have a 
greater discriminatory power in isolating the NB and FB 
phases than X-ray fluxes alone. However, even optical fluxes 
alone cannot properly separate the two phases. 

(ii) The PSD of Sco X-l appears to change in shape be¬ 
tween low and high optical luminosity states. Qualitatively, 
the low-luminosity PSD can be expressed as a red noise PSD 
with a power-law ~ —2. During high optical luminosities 
an additional component is observed on top of the power- 
law with a possible break at ~ 10 -3 Hz. This additional 
broad-band feature could possibly be associated to unstable 
nuclear burning on the NS surface (Bildsten 1993). 

(iii) Sco X-l displays rms-flux relations in both the low 
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and high optical luminosity states. However a linear rela¬ 
tion is only observed during low luminosities. The deviation 
from linearity, during high luminosity states, might suggest 
that an additional variability component is present, thus re¬ 
moving the linearity and stationarity at the highest fluxes. 
This component could possibly be associated with the extra 
variability component observed in the FB PSD. 

(iv) We find that the equivalent widths of the low ioni¬ 
sation lines Ha 6563 A, H f3 4861 A, Hel 5876 A, Hel 6678 
A and Hel 7065 A are higher during times of lower optical 
luminosity, whilst Hell 4685 A does not appear to change 
in strength. This can be explained if during high optical 
luminosity periods (when M is highest) the accretion disc 
is too hot for the production of low ionisation lines, but 
not hot enough to remove higher ionisation lines such as 
Hell 4685 A. We also find that Hell 4685 A appears to be 
double-peaked only during high optical luminosity states. It 
is possible that the single-peaked profile of Hell 4685 A is 
associated to an accretion disk outflow during the low op¬ 
tical luminosity phase. This interpretation is in line with 
previous X-ray measurements of the Iron K a emission line 
model (Bradshaw et al.| 2007), where fast outflows are sup¬ 
posed to be generated in Sco X-l during periods of high 
mass transfer. 

(v) The DCF between the available optical K2 and X-ray 
MAXI lightcurves displays clear asymmetry. The DCF pro¬ 
file displays a clear, broad, 12.5 hour secondary optical lag 
peak. This can potentially be attributed to the disk thermal 
reprocessing timescale. 

(vi) The DCF obtained during the FB also displays clear 
asymmetry, and is qualitatively similar to the full DCF, with 
a secondary marginal detection of a 4.5 hour optical lag 
peak. This too could potentially be the signature of the disk 
thermal timescale. It is thus possible that the broad optical 
lags are mostly generated during the FB state. 

(vii) The DCF obtained during the NB clearly show that 
the optical and X-ray variability are anti-correlated during 
this phase. This phenomenon has been previously reported 
by other authors. 

(viii) We detect clear echoes within the optical K2 light 
curve alone using the Cepstrum, at 1 hour, 4.4 hours and 
12.6 hours. The 12.6 hour echo timescale is consistent with 
the full optical/X-ray DCF, whilst the 4.4 hour echo is con¬ 
sistent with the FB DCF. Even if it were present in any of 
the DCFs studied here, we would not be able to recover the 
1 hour echo due to our DCF resolution. We interpret these 
echoes as the disk thermal timescale, connected to the DCF 
optical lags. Thus, as some of the X-ray radiation will be re¬ 
flected quasi-instantaneously by the outer optical emitting 
disk, some radiation will be absorbed and re-emitted on the 
thermal timescale. This would give rise to the broad opti¬ 
cal lags observed. The reprocessing mechanism would also 
generate echoes within the optical lightcurve alone, with a 
delay close to the thermal timescale. 
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